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In  this  paper,  the  modelling  of  an  energy  generation  system  based  on  polymer  electrolyte  membrane  fuel 
cell  (PEMFC)  system  through  a  parameter  varying  approach  (LPV  model),  that  takes  in  to  account  model 
parameter  variation  with  the  operating  point,  is  presented.  This  model  has  been  obtained  through  a 
Jacobian  linearization  of  the  PEMFC  non-linear  dynamic  model  that  was  previously  calibrated  using  real 
data  from  lab.  In  order  to  illustrate  the  use  of  the  LPV  model  obtained  its  application  to  model-based  fault 
detection  is  used.  For  this  purposes  a  set  of  common  fault  scenarios,  which  could  appear  during  a  normal 
PEMFC  operation,  is  used  as  case  study. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Low-temperature  PEM  fuel  cells  are  considered  as  sources  for 
rapid  medium  of  energy  generation,  making  these  equipment  suit¬ 
able  for  automobile  applications.  The  supply  of  raw  materials 
(usually  air  or  pure  oxygen)  is  normally  performed  using  an  air 
compressor  or  blower  and  hydrogen  stored  in  tanks.  The  system 
uses  additional  equipment  to  carry  materials  reaction  to  the  opti¬ 
mum  operating  conditions,  such  as  cooling  systemsand  humidifier. 
During  the  chemical  reaction  that  is  taking  place  into  the  stack, 
where  the  energy  is  generated,  different  phenomena  occur,  such  as 
thermal,  fluid-mechanical  and  electrolytic. 

The  complex  and  non  linear  dynamics  of  the  power  generation 
systems  based  on  fuel  cell  technology  lead  to  the  use  of  linear  mod¬ 
els  that  includes  parameter  varying  with  operating  point  not  only 
for  advanced  control  techniques  but  also  for  fault  diagnosis  algo¬ 
rithms  based  on  models.  The  use  of  LPV  models  is  an  alternative 
to  the  approaches  presented  in  previous  works  [1,2]  addressing 
methodologies  for  monitoring  and  fault  diagnosis  based  on  a  the¬ 
oretical  non-linear  dynamic  model  proposed  by  Pukrushpan  [3,4]. 
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Within  the  recent  decade,  state  of  art  and  background  about 
control  of  LPV  systems  has  been  developed  [5-10].  Because  of  a 
LPV  system  can  be  considered  as  a  parametrized  family  of  linear 
systems  that  change  with  the  operating  point  conditions,  then  LPV 
technique  allows  a  systematic  approach  for  control  and  fault  diag¬ 
nosis  system  design.  At  the  cost  of  conservatism  the  approach  can 
be  applied  to  an  even  wider  range  of  systems  known  as  quasi- 
LPV  systems,  where  varying  parameters  are  scheduled  with  state 
variables. 

Since  LPV  models  are  structured  as  similar  as  a  linear  time- 
invariant  (LTI)  state  space  system,  the  control  and  fault  diagnosis 
design  methods  can  easily  be  extended.  The  main  contribution  of 
this  paper  is  to  obtain  a  linear  parameter  varying  model  for  a  typical 
PEMFC  and  illustrate  its  use  for  robust  fault  detection  using  interval 
observers. 

2.  Fuel  cell  modelling 

The  model  proposed  in  [1 1  ],  is  a  non-linear  dynamic  model  cali¬ 
brated  using  real  data  from  laboratory  using  a  Isq- non  linear  fitting 
approach  [12,13].  This  model  is  able  to  reproduce  the  behaviour 
of  a  commercial  PEMFC  (Ballard  1.2  kW,  Nexa®)  prototype,  which 
has  been  identified  in  a  wide  range  of  operating  conditions.  Fig.  1 
shows  the  dynamic  model  layout. 

2.1.  Dynamic  non-linear  model 

The  model  is  considered  as  SIMO  system,  where  the  input  (u) 
is  the  stack  current  (lst)  and  the  outputs  (y)  are  battery  tempera¬ 
ture  (Tst),  stack  voltage  ( vst ),  oxygen  consumption  ratio  (Aq2  ),  speed 
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Fig.  1.  Nexa®  PEMFC  simulator  schematic. 


engine  (cocp)  and  inlet  pressure  to  the  cathode  (Psm).  The  voltage 
from  the  compressor  (vcm)  is  controlled  using  a  static  feed-forward 
controller.  Fig.  2  shows  the  layout  of  each  subsystem  embedded 
into  the  PEMFC  dynamic  model. 

The  proposed  model  consists  often  state,  and  the  state  equations 
are  listed  in  the  following 

&>cp  =  i - ('Lem  ~  Pep), 

JcpLOcp 

Prm  =  (Wra, o  -  Wrm, o), 

Vrm 

rhrm  =  Wca  ~  Wrm,o, 

Psm  =  WcpTcp  -  Wsm,0Tsm), 

Vsm  (1) 

131  sm  =  WCp  —  VVsm,o? 

™H2  =  WH2,i  -  ^H2,o  -  WH2,r-WH2,nr 
™w,an  =  WVan,i  ~  WVan,o  ~  WVmbr , 

1TIN2  =  WNl j  -  Wn2,0, 

irio2  =W0 2J  -  Wo2,0  -  Wo2,r, 

IPlstCstTst  =  Hreac  ~  P elec  ~  Qrad  ~  Qconv- 

The  state  variables  (x)  of  this  dynamic  model  are  the  following: 
mass  of  oxygen  (m02),  nitrogen  (mN2),  hydrogen  (mH2),  cathode 
water  flow  (mw?Ca)»  stack  temperature  (Tst),  angular  velocity  of  the 
compressor  ( ooCp ),  supply  pressure  ( Psm )  and  return  pressure  (Prm) 
of  the  humidifier,  inlet  flow  ( msm )  and  outlet  flow  (rrirm)  of  humid¬ 
ifier.  The  subindex  in  the  variables  i,  o,  r,  ml  means,  input,  output, 
reaction  and  natural,  respectively.  In  the  heat  balance  the  subindex 
reac ,  elec ,  rad  and  conv  are  related  respectively  to  reaction,  electric, 
radiation  and  convection. 

The  system  perturbation  (z)  that  have  been  considered  are 
related  to  the  weather  conditions  ( Tamb ,  Patm\ 

The  model  output  equations  are: 


•  Compressor  speed  motor: 


where  vst,  total  stack  voltage  (V);  E,  open  circuit  voltage  (V); 
Pact,  activation  voltage  loss  (V);  voh ,  ohmic  voltage  lose  (V);  pcon, 
concentration  voltage  loss  (V);  rifC ,  amount  of  cells;  A,o2,  oxygen 
excess  ratio;  /st,  stack  current  (A);  F,  Faraday  constant  (col/mol); 
W,  mass  flow  (g/s);  Ucp,  compressor  blade  (KRPM);  dc,  compres¬ 
sor  diameter  blade  (%);  p,  compressor  efficiency  (%);  Mq2  ,  oxygen 
molar  weight  (g/mol);  Mh2  ,  hydrogen  molar  weight  (g/mol);  x02 , 
oxygen  fraction  (%);  (pit  humidity  (i  =  ca,  an)  (%). 

3.  Linear  parameter  varying  model 

Exist  different  ways  to  obtain  LPV  models.  Some  methods  use 
non-linear  equations  of  the  system  to  derive  a  LPV  model  such  as 
state  transformation,  substitution  of  functions  and  methods  using 
the  well  known  Jacobian  linearization  [14-16].  Another  kind  of 
method  uses  multi-model  identification  that  consists  basically  in 
two  different  steps:  (1)  a  set  of  LTI  model  is  identified  at  differ¬ 
ent  equilibrium  points  by  classical  methods  ( on-line  or  off-line),  (2) 
then,  the  following  step  is  to  get  a  multi-model  by  an  interpola¬ 
tion  law  that  allows  to  commute  among  local  LTI  model  at  each 
operating  point  [17,18]. 

3.1.  Problem  formulation 

The  type  of  LPV  system,  which  is  considered  in  this 
paper,  assumes  an  affine  dependence  with  a  parame¬ 
ter  vector  &k  and  can  be  described  in  discrete  time  state 
space  as: 


•  Stack  voltage: 

Pst  =  Life  '  (£  —  Pact  —  Poh  —  Peon)- 

•  Oxygen  excess  ratio: 

.  _  Wo 2,i  _  *o2  •  Wcp 

°2  “  Wo^  -  (Mo2  n/c./st)/4.F- 


xk+1  =  A(tffc)xfc  +B(^)uk  +  wfc, 

(2)  Yk  =  C(4)xk  +  D(&k)uk  +  vk, 

where  xkeRnx,  U/{Glnu  and  y keRny  are,  respectively,  the  state, 
input,  and  output  vectors.  The  process  and  measurement  noises  are 

(3)  eRn*  and  vk  e Rny  respectively.  Both  are  considered  unknown 
but  bounded  as  vk  e  Yn x  and  e  W  which  are  interval  boxes. 
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Fig.  2.  Nexa®  PEM  fuel  cell  dynamic  model  diagram. 


&ke®  is  the  vector  of  time-varying  parameters  that  change 
with  the  operating  point  scheduled  by  some  measured  system 


function: 

h  =f(Pk), 
where 
0  =  {&ke 


<h<  h)- 


The  system  described  in  Eq.  (5)  could  be  seen  as  a 


tope  of  dimension  is  N  described  as  an  array  of  matrices 
described  by 


Hh) 

CW 


B(4) 

Wk) 


Dj(i9J ) J  J 


j=i 


B j(&) 

V0) 

A 0) 
C ){&) 


Bj(^) 

Dj(^) 


(8) 


where  Ay,  By 


Cyand  Dj 


are  the  state  matrices  for  each  jth  model 
obtained  by  linearization  of  the  non-linear  model  around  the  jth 
operating  point.  Using  Eq.  (8),  the  Eq.  (5)  can  be  approximated  as 
follows: 


Xfc+1  =  ^V(pk)  •  [Aj(#0  B,(!?J)]  • 


i=l 

N 


U  k 


(9) 


yk  =  ^(pk)  •[(*(#)  D0)]- 

j= i 

3.2.  PEMFC  LPV  structure 


Using  the  PEMFC  non-linear  dynamic  model  presented  in 
Section  2.1  and  using  the  Jacobian  linearization  approach,  the  non¬ 
linear  model  is  possible  to  transform  it  into  a  LPV  model  in  state 


space  form,  as  in  Eq.  (9)  with  the  following  system  matrices  struc¬ 
ture 


"ail 

0 

0 

014 

0 

0 

0 

0 

0 

0 

0 

«22 

«23 

a24 

«25 

0 

0 

«28 

«29 

0210 

(6) 

0 

«32 

0 

a34 

035 

0 

0 

a38 

039 

a310 

a4i 

0 

0 

a44 

«45 

0 

0 

«48 

a49 

a4i0 

A,  = 

a5i 

0 

0 

054 

0 

0 

0 

ass 

059 

asio 

J 

(7) 

0 

0 

0 

a64 

0 

a66 

a67 

0 

0 

a6io 

0 

0 

0 

a74 

0 

a76 

a77 

0 

0 

0710 

lin- 

0 

a82 

0 

«84 

ass 

0 

0 

ass 

«89 

a810 

Eq. 

0 

«92 

0 

a94 

095 

0 

0 

ag8 

agg 

agio 

olv- 

_  0 

0 

0 

0 

0 

a106 

aio7 

0 

a109 

0 

C,  = 


Bj  = 


Cn  0  0 

0  0  0 

0  0  0 

0  0  0 


b ii  b 12 
0  0 


0 

C24 


c34  C35 


0  0  C45  C47 


C38  c39  C310 

C48  C49  C410 


0 

0 

0 

bn 

^81 

0 

^91 

^101 


DJ  = 


0  0 
0  0 
d3i  0 
d4i  0 


where 


(10) 


•  States:  X-  [C0cp  Pom  Worn  Psm  ITlsm  a?H2  ^w,an  ^02  Tst]T. 

•  Inputs:  u  =  [/st  Vcm]T  where  the  scheduling  variable  is  Jst- 

•  Outputs:  y  =  [Psm  vst  A q2  cocp]. 

•  Perturbations:  z=[Tamb  Patm\. 

In  this  case  the  scheduling  variable  is  the  current  (Jst).  Thus  Eq. 
(6)  is  =f(Istk )  and  the  varying  parameters  (ay,  by,  Cy  and  dy) 
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of  the  system  matrices.  Consequently  the  system  in  Eq.  (9)  can  be 
particularized  as  follows. 

N 

Xfc+1  =  y>(/sg)  ■  [Aj(ffl) 
i= i 

N 

yk  =J20liVst,k)-[Cjm 

J=1 

4.  Linear  parameter  varying  observer 

4.1.  Definition 

Using  the  PEMFC  LPV  model  obtained  in  previous  section,  a 
discrete  time  LPV  observer  for  state  estimation  with  Luenberger 
observer  structure  can  be  expressed  as 

Xfc+1  =  A  (&k)xk  +  B(&k)uk  +  L(4Xyk  -  yk)  +  wk,  (u) 

y k  =  C(#k)xk  +  D(£k)uk  +  vk. 

Using  Eq.  (1 1 )  the  observer  can  be  alternatively  expressed  as 

N 

Xk+1  =  ^(/sa)[A“(^)  . 

j= 1 

N 

yk  =  p <,(&)]  ■ 

J=1 


*k 

*k 


(11) 


where 

A-  =  Aj(^)  -  Lj(^)  •  Cj(^). 

The  observer  gain  (L)  should  be  designed  to  guarantee  closed  loop 
stability  for  all  operating  points,  i.e.  &k  e  S.  This  is  achieved  through 
LMI  formulation  for  pole-placement  within  a  wide  class  of  pole 
clustering  regions,  defined  by  an  affix  {-q,  0)  and  a  radius  r  such 
that  (-q  +  r)<0  [22]. 


4.2.  Observability  analysis 


The  observer  design  implies  the  fulfilment  of  the  observability 
condition  that  is  related  to  the  number  of  states  that  can  be  inferred 
(or  estimated)  based  on  the  available  measured  outputs.  For  the 
study  of  the  system  observability  degree  at  each  operating  point 
(#)  the  following  conditions  are  used 


rank 


Xfi  -  A 

C 


cond 


(M)- 


(14) 


As  discussed  in  [4],  a  given  system  eigenvalue  is  unobservable 
if  the  rank  condition  indicates  a  rank  loss.  On  the  other  hand,  a 
large  condition  number  implies  that  the  corresponding  eigenvalue 
is  weakly  observable. 


5.  Application  to  fault  detection 

In  order  to  test  the  proposed  modelling  methodology  for  fault 
detection,  a  common  set  of  faults  have  been  selected  as  benchmark. 
This  set  of  fault  was  included  in  the  model  (Ballard,  Nexa®)  [1 1  ]  in 
simulation  environment  in  the  MATLAB/SIMULINK®  (see  Fig.  1). 

5.1.  Model  analysis 

The  stability  of  the  PEMFC  system,  described  in  LPV  form  by 
Eq.  (5),  at  each  operating  point  is  verified  if  all  the  eigenvalues  of 
the  system  are  negative.  In  the  case  of  the  a  operating  point  corre¬ 
sponding  to  1st  =  15  A,  the  stability  of  the  system  can  be  verified  by 


Table  1 

Eigenvalues,  rank  and  condition  number  for  the  observability  analysis  of  the  poly¬ 
topic  LPV  system  at  Ist  =  1 5  A. 


State 

A., 

Rank 

Cond 

OJcp 

-1.00E+04 

10 

1.11E-03 

Prm 

-543.8 

10 

27.92 

mrm 

-388 

10 

31.87 

Psm 

-49.9 

10 

404.82 

msm 

-13.14 

10 

2.22E+03 

mH  2 

-4.88 

10 

1.04E+04 

TO  an 

-1.72 

10 

5.23E+04 

mNz 

-1.26 

10 

7.53E+04 

mo2 

-0.003 

10 

1.90E+07 

Tst 

-2.50E-07 

10 

3.31E+11 

Table  2 

Eigenvalues,  rank  and  condition  number  for  the  observability  analysis  of  the  reduced 
polytopic  LPV  system  at  Ist=  15  A. 

State 

A.f 

Rank 

Cond 

-26.81 

4 

30.4 

*2 

-4.91 

4 

11.5 

*3 

-2.21 

4 

88.94 

*4 

-3.50E-03 

4 

3.80E+03 

the  eigenvalues  shown  in  Table  1. 

On  the  other  hand,  although  the 

system  is  observable  using  the  observability  matrix  test,  the  observ¬ 
ability  of  the  system  could  be  further  analyzed  with  conditions 
showed  in  Eq.  (14). 

From  Table  1,  it  is  clear  to  see  that  there  are  states  ( cocp ,  Psm, 
irism ,  Prm,  knrm )  which  contain  fast  dynamics  compared  to  others 
such  as  the  temperature  (Tst)  or  oxygen  mass  (m02).  Both  types 
of  dynamics  could  decrease  the  observability  performance  of  the 
system  producing  problems  with  the  state  estimation  using  the 
observer  (Table  2).  This  can  be  verified  by  looking  at  condition 
number  introduced  in  Eq.  (14).  To  solve  this  problem,  the  order 
of  the  model  should  be  reduced  as  proposed  in  [23].  Table  1  shows 
that  the  eigenvalues  of  the  four  first  states  present  a  higher  value 
than  the  rest  of  them.  This  means  that  the  model  can  be  reduced  to 
a  four  order  state  space  model  instead  of  the  original  ten  order 
state  space.  The  reduced  system  is  observable  and  the  condi¬ 
tion  number  has  improved  with  respect  to  the  values  presented 
in  Table  3. 

In  order  to  test  the  behaviour  of  the  reduced  model  compared  to 
the  complete  model,  both  models  were  implemented  in  simulation. 
Fig.  3  shows  the  simulation  results  of  the  complete  and  reduced 
PEMFC  LPV  model  under  a  sequents  of  operating  point  changes  in 
current  demand  (Jst).  This  figure  shows  also  the  measured  outputs 
for  comparison  with  those  generated  in  simulation.  Fig.  4  shows 


Table  3 

Description  of  the  fault  scenarios  implemented  in  FGB. 


ID 

Fault  description 

Type 

Magnitude 

/i 

There  is  a  suddenly  increase  of 
friction  in  the  mechanical 
component  part  of  the  compressor 

Parametric  abrupt 

AkR  =  60% 

h 

Degradation  in  the  stack  cells  stack 
is  presented  because  of 
contact-sensitivity  reactions 
against  to  a  reaction  killer 

Parametric  abrupt 

40% 

fs 

Hydrogen  leak  in  the  anode  is 
presented  because  of  seal 
degradation. 

Parametric  abrupt 

An//= 2E+2 

/4 

A  suddenly  leak  of  hydrogen  is 
presented  at  the  anode  inlet 
manifold 

Parametric  abrupt 

80% 

fs 

A  suddenly  leak  of  air  is  presented 
at  the  inlet  outlet  supply  manifold 
inlet  manifold 

Parametric  abrupt 

10% 
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Fig.  3.  Comparison  of  process  output  signal  between  real  and  reduced  polytope  LPV  models. 


the  model  error  between  the  models  and  the  measured  outputs 
of  the  real  process.  From  either  of  the  two  previous  figures,  it  can 
be  noticed  that  the  reduced  model  presents  better  performance 
than  the  complete  one.  This  result  justifies  the  model  reduc¬ 
tion  since  improves  the  observability  and  numerical  conditioning 
properties. 

5.2.  Fault  benchmark 

The  usefulness  of  LPV  model  and  observer  proposed  in  this  paper 
is  illustrated  in  fault  detection.  For  this  purpose,  a  set  of  common 
possible  fault  scenarios  is  considered  and  implemented  using  the 
simulator  presented  in  Fig.  1  by  adding  a  fault  generator  block  (FGB) 
( see  Fig.  7).  Table  3  describes  the  set  of  faults  which  were  considered 
as  case  study.  In  the  following  section,  it  is  described  how  the  faults 
were  included  in  simulation. 


5.2.1.  Fault  1 

The  fault /i  is  simulated  with  an  increment,  A Rcm,  in  the  com¬ 
pressor  motor  resistance  Rcm •  The  fault  effect  is  translated  in  a 
change  in  the  compressor  torque  rcm  affecting  directly  the  state 
variable  aoCp 


1cm 


Hcmkt  (ttt 
Rcm  +  Ai?cm  Cm 


kv(L>cm), 


(15) 


where  r\cm  is  the  motor  mechanical  efficiency,  kt  and  kv  are  motor 
constants,  and  ooCm  is  the  compressor  speed.  Furthermore,  the 
parameter  is  related  with  the  state  of  Psm  because  its  dynamics  are 
governed  by  the  compressor  inlet  air  flow 


Wca,i  —  (ksm  +  A kSm){Psm  ~  Pep),  (16) 

where  Pcp  is  the  compressor  pressure  and  Psm  is  the  supply  manifold 
pressure.  The  amount  of  flow  that  is  feed  into  the  cathode  is  related 


x  10" 
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Error  en  r 


Error  en  V, 


St 


02 


Error  en  oo 


cp 
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Fig.  4.  Model  residuals  (reduced  model  and  real  process  data). 
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- Upper  Boundary  - Lower  Boundary  -  Process  output 


Fig.  5.  Confidence  intervals  for  predicted  outputs  against  real  measurements  in  faultily  operation  (/2,/4). 


to  pressure  with  a  linear  constant  ksm-  It  is  possible  to  simulate  a 
reduction  in  the  supply  manifold  as  a  change  in  ksm  as  Aksm-  Thus, 
the  air  mass  flow  is  affected  by  the  fault  and  consequently,  the  total 
mass  balance  across  the  PEMFC  changes. 


as  Anijr-Anlf3. 

A nU 

'\/  Ranlpa 


wH2,n;  = 


)V2 


(19) 


5.2.2.  Fault  2 

Fault  f2  is  presented  as  any  contamination  in  the  reaction  into 
the  stack  reducing  the  chemical  reaction  efficiency  by  decreasing 
of  catalysis  active  area.  The  current  density,  i,  is  defined  as  current 
per  area  in  a  single  cell,  which  is  equals  to  stack  current  Jst  (A),  per 
cell  active  area,  AfC  (cm2).  If  the  stack  is  mead  of  a  set  of  cells  in 
series  architecture. 


i  = 


ht 


Afc  •  AA/c  ’ 


(17) 


where  AAfC  is  the  active  contaminated  area. 

Because  of  the  major  voltage  drops  are  associated  with  cur¬ 
rent  density  for  non-linear  relations,  see  [24],  current  density  is 
an  important  issue  for  total  stack  voltage 


where  Pr  =  PanlPatm  is  the  pressure  ratio  across  the  assumed  leak 
and  the  anode  gas  constant,  Ran,  is  calculated  through  the  universal 
gas  constant,  R  as  follows 

an  yH2MH2  +(1  -yH2)MH2o’  ^  ) 

where  the  molar  fraction  of  hydrogen  in  the  anode  is  given  by 


y«2 


Pan 


<PanPsatj=st 

Pan 


(21) 


5.2.4.  Fault  4 

This  fault  introduces  a  leak  of  hydrogen  and  is  simulated  as  a 
change  in  the  mass  balance  in  the  hydrogen  inlet  flow  as  follows 


Pst  —  Life  x  [F  —  Pact  ~  Pohm  ~  Vconc],  (18) 

where  E  is  the  open  circuit  voltage;  vact,  v0hm  and  vConc  are  activation 
loss,  ohmic  loss  and  concentration  loss,  respectively,  then  is  clear 
that  fault/2  will  have  a  direct  action  over  the  output  of  stack  voltage. 

5.2.3.  Fault  3 

The  term  nt  introduced  in  Eq.  (10)  represents  the  natu¬ 
ral  leak  from  the  anode  of  the  fuel  cell  stack.  This  leak  is  always 
present  due  to  the  physical  stack  sealing  design.  It  is  assumed  that 
the  natural  leak  is  governed  by  a  standard  orifice  relation  through 
an  effective  area,  Ani.  This  parameter  has  been  obtained  in  [25].  In 
order  to  simulate  a  degradation  in  the  seal  a  change  in  An{  is  used 


wH2,f  =  wH2,i-WH2,if4.  (22) 

5.2.5.  Fault  5 

This  fault  appears  as  a  leak  of  air  at  the  cathode  inlet  flow. 
Because  this  fault  is  considered  as  a  leak,  it  is  introduced  in  the 
mass  balance,  as  fault /4,  as  follows 

WhmJ  =  WhmJ-WhmJf5.  (23) 

Note  that  the  amount  of  air  that  does  not  enter  into  the  system  will 
not  only  create  an  abrupt  change  in  the  total  mass  balance,  where 
mo2 ,  mN2  state  variable  are  mainly  involved,  but  also  a  system  pres¬ 
sure  change. 
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Fig.  6.  Confidence  intervals  for  predicted  outputs  against  real  measurements  in  faultily  operation  C/2,/4). 


5.3.  Fault  detection  systems 

The  fault  diagnosis  system  (FDS)  based  on  the  LPV  observer 
proposed  in  this  paper  has  been  implemented  using  MAT- 
LAB/SIMULINK  environment  [26].  The  simulator  of  fuel  cell  model 
(Ballard  Nexa® )  presented  in  Section  2.1  is  used  as  virtual  reality  for 
the  fault  detection  case  study.  The  FGB  and  FDS  subsystems  blocks 
were  added  to  the  PEM  FC  simulator  to  create  the  fault  scenarios 
presented  in  the  previous  section. 

5.3.1.  Fault  detection  process  (FDP) 

Using  the  measured  inputs  and  outputs  presented  in  Section  2.1 
and  using  the  structural  analysis  methodology  [27],  the  following 
set  of  residuals  can  be  obtained  as: 

n  =  Psm  -  Psm , 

r2  =(t>cp-d)Cp,  (24) 

*3  =  Ao2  -  A,q2  , 
r4  =  Vst  -  Vst  • 


The  residuals  are  defined  conceptually  as  the  differences  the  pro¬ 
cess  measurements  and  output  estimation  from  the  LPV  observer 
as  described  in  Section  3.2.  From  the  model  analysis  in  Section  5.1, 


the  reduced  model  offers  a  better  performance  when  implement¬ 
ing  the  LPV  observer  than  the  complete  model.  Thus  the  reduced 
model  is  used  for  creating  the  PEMFC  LPV  model. 

Fig.  5  shows  the  measured  and  the  confidence  intervals  for  pre¬ 
dicted  outputs  considering  noise  and  process  noise.  Notice  that 
when  the  /2  is  introduced  at  k  =  50  in  the  PEMFC  system,  at  the 
time  /<  =  85  the  sensor  measurements  cross  the  boundaries  of  con¬ 
fidence  intervals  of  predicted  outputs.  This  allows  detecting  the 
fault. 

Fig.  6  shows  PEMFC  process  behaviour  compared  with  the  adap¬ 
tive  thresholding  generated  the  interval  observer.  The  process 
suffers  a  degradation  in  supply  pressure  (Psm)  sensor  fault  over  the 
time  starting  from  k  =  20  and  ends  in  k  =  60,  the  measurement  cross 
the  upper  boundary  at  /<  =  40.  Note  that  the  other  variables  do  not 
cross  their  thresholds.  Later  appears  a  sensor  offset  in  Xq2  ,  that  leads 
the  measurement  to  cross  the  upper  bound  at  /<  =  60.  The  effect  of 
this  second  fault  additionally  acts  over  the  interval  estimation  of 
®cp  (Fig.  7). 

The  confidence  intervals  for  predicted  outputs  are  computed 
using  zonotope-based  algorithm  proposed  in  [28],  which  offers  an 
efficient  way  of  taking  into  account  model  uncertainty  when  esti¬ 
mating  the  state  using  an  LPV  observer. 


Nexa  PEMFC  non  linear  model 


xk+i  =  A(d)xk  +  B(ti)uk  +  --- 
=  L(d)  (yjt-yjk) 
yjfc  =  C  (fl)  Xjt  +  D  (&)  uk 


y/ 


Fig.  7.  Fault  generator  block  (FGB)  and  fault  diagnosis  system  (FDS)  implementation  diagram. 
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6.  Conclusions 

In  this  paper,  an  LPV  model  is  introduced  as  a  means  to  approx¬ 
imate  a  PEMFC  that  presents  a  highly  non-linear  behaviour  using 
a  model  with  linear  structure  but  with  parameters  that  vary  with 
the  operating  point.  This  model  can  be  obtained  from  the  lineariza¬ 
tion  of  non-linear  model  which  has  been  calibrated  above  a  set  of 
data  from  a  commercial  PEMFC  model  (Ballard,  Nexa®)  around  a 
set  of  operating  points.  In  addition  to  the  LPV  model,  the  paper 
proposes  the  design  of  an  LPV  observer  to  estimate  PEMFC  system 
states.  To  illustrate  the  usefulness  of  the  LPV  model  and  observer, 
the  application  to  fault  detection  is  used.  The  paper  also  analysis 
the  observability  of  the  LPV  model  obtained  concluding  that  the 
complete  LPV  model  contains  dynamics  that  are  very  different  that 
create  numerical  difficulties  in  the  observation  implementation. 
These  difficulties  can  be  overcome  by  model  reduction  neglecting 
system  fast  dynamics.  The  application  of  the  LPV  observer  to  fault 
detection  is  shown  using  a  simulator  developed  using  an  exper¬ 
imentally  calibrated  nonlinear  model  of  a  typical  PEMFC  (Nexa® 
Power  Module).  A  set  of  common  fault  scenarios  have  been  defined 
and  implemented  in  the  simulator  to  test  the  LPV  observer  in  fault 
detection.  Finally,  the  paper  presents  how  the  observer  can  satis¬ 
factorily  detect  faults  in  some  of  the  fault  scenarios  defined.  As  an 
extension  of  the  research  presented  in  this  paper,  an  algorithm  that 
allows  not  only  detect  faults  but  also  isolate  them  is  being  devel¬ 
oped  based  also  in  the  LPV  observer  developed.  Some  preliminary 
results  have  already  been  obtained  that  show  promising  results  in 
this  line. 
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